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ABSTRACT 

Using  a  rotating  cone,  a  study  was  made  to  determine  the  velocity 
dependent  shear  property  of  various  soils.  In  the  wide  variety  of  soi 's 
tested  Including  clays,  wet  and  dry  sands,  no  discernible  velocity  dependent 
shear  properties  were  measured 
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SOIL  STUDIES  -  AN  INVESTIGATION  INTO  THE  PROPERTIES  OF  THE  SOIL-WHEEL  INTERFACE 


Part  2  -  Results  of  Tests  with  a  Rotating  Cone  in  Sand 

by 

Lou  is  I •  Le v i t i cus 

INTRODUCTION 

The  purpose  of  this  investigation  was  to  verify  the  existence  of 
velocity  dependent  soil-wheel  interface  parameters  as  postulated  and 
developed  in  the  first  part  of  this  study. '  The  basic  premises  for  the 
theoretical  considerations  were: 

a.  There  exist  soil-rigid  bodies  interface  properties  which 
manifest  themselves  when  relative  moment  occurs. 

b-  The  interface  properties  arc  velocity  dependent- 

c-  The  interface  properties  can  be  separated  into  cohesive 
and  frictional  properties. 

In  order  to  determine  these  properties,  a  mechanism  is  needed  which 
can  generate  relative  movement  of  the  rigid  body  with  respect  to  the 
soil  mass  at:  various  velocities. 

Moreover,  it  is  required  that  the  phenomenon  as  it  occurs,  will  have 
sufficient  similarity  to  the  phenomena  occurring  near  a  moving  soil-wheel 
interface  to  warrant  the  use  of  those  parameters. 

Due  to  financial  restraints,  the  testing  was  performed  mainly  on  sand. 
The  tests  performed  on  clay  were  not  numerous  enough  to  resul  :  in  reliable 
conclusive  data  and  are  therefore  not  included  in  this  report- 
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During  the  execution  of  this  work  a  study  was  published  on  soil 

2 

wheel  interaction  at  high  speeds  •  The  author  conducted  tests  of  full 
size  aircraft  landing  gear  wheels  at  velocities  up  to  about  100  knots. 

The  soils  used  were  a  buckshot  clay  and  a  rivet  sand.  The  tests  showed 
a  reduction  of  drag  load  and  rut  depth  with  increasing  velocity  up  to  about 
20  knots-  After  this  minimum  the  drag  load  and  rut  depth  increased  to  a 
maximum  in  the  vicinity  of  kO  knots-  It  was  tentatively  assumed  that 
up  to  20  knots  no  inertial  effects  took  place.  Below  this  threshold 
one  would  expect  that  the  viscous  effects  dominate.  The  tests  conducted 
with  the  cone  were  at  lower  velocities,  up  to  about  10  knots-  Also  our 
tests  were  so  far  conducted  only  cn  a  sand-  Thus  our  data  are  not  directly 
comparable  with  the  data  in  the  paper. 


The  author  used  a  modified  version  of  the  WES  sand  and  clay  mobility 
numbers-  The  modification  does  result  in  two  dimensional  mobility 


numbe  rs  • 


0 

c 


Cl  (bd) 


where  Cl  ~  Cone  Index  (psi) 

G  =  Cone  index  gradient  (psi/in) 
l  =  Wheel  width  (in) 
d  =  Wheel  diameter  (in) 

=  Drag  force  (lb) 
h  =  Tire  section  height  (in) 
tj  =  Tire  deflection  (in) 
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These  dimensions  may  be  part  of  the  reason  for  some  of  the  differences 
between  the  theoretical  and  empirical  curves  since  'J.  and  are 

both  used  in  the  rut  depth  prediction  equation  and  in  the  drag  force 
equa  t ion . 

The  nature  of  the  empirical  drag  force  curves  in  the  first  region 
for  a  sandy  soil  do  show  a  similar  trend  as  did  the  torque  versus 
velocity  data  for  the  rotating  cone. 

3 

Additional  work  was  performed  at  WES  by  Turnage  who  penetrated 
cones  at  different  velocities  in  the  soil.  It  is  not  clear  al  this 
point  how  the  results  of  our  rotation  tests  can  be  compared  to  the 
penetration  tests. 
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Design  of  the  Instrument 

Tlie  design  consisted  of  a  cone  mounted  on  a  shaft  and  driven  by  a 
servo-controlled  motor  (Figs-  I  ,2,3)  * 

Drive  motor,  bearings,  cone  and  velocity  transducer  are  mounted  on 
a  plate,  which  is  attached  to  a  second  plate  by  means  of  a  set  of  radial 
arms  which  allow  small  radial  movement-  This  movement  is  sensed  by  a 
Shaevi tz  LVDT  and  calibrated  with  a  known  torque  input  (Figure  2). 

The  force  of  penetration  is  measured  by  a  Lebow  load  cell  mounted 
below  the  sample  tray-  A  position  crank  is  used  to  zero  the  sinkage 
measuring  potentiometer  and  the  sample  is  pulled  upward  along  the  shafts 
against  the  cone  by  the  weights  in  the  two  weight  pans-  This  type  of 
operation  was  developed  because: 

a-  The  cone  finds  its  own  depth  under  the  actual  applied  weight- 

b.  The  accuracy  of  measurement  is  improved  since  the  weight 
of  the  soil  sample  is  much  less  than  that  of  the  cone,  drive  motor  and 
torque  transducing  system- 

c-  The  load  on  the  sample  can  be  recorded  continuously 
while  the  tes  t  is  run . 

d-  The  removal  of  the  sample  and  insertion  of  a  new  sample 
is  made  easier- 

e-  The  drive  system  and  transducers  could  be  mounted  more 
rigidly  on  the  framework  of  the  machine. 

Two  stainless  steel  cones  were  made  having  apex  angles  of  30° 
and  60°.  The  30°  cone  was  5"  high  and  the  60°  cone  was  4"  high. 
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FLOATING  PLATE- 
CARRIES  CONE  AND' 
DRIVE 


FIG.  2.  DETAIL  OF  TORQUE  BALANCE 
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The  following  quantities  were  measured  during  a  test' 

1.  the  penetration  force  on  a  cone 

2.  the  torque  needed  to  rotate  the  cone 

3.  the  rotational  velocity 

4.  the  sinkage  of  the  cone 

Preliminary  calculations  indicated  the  following  ranges: 

1.  torque  values  up  to  6  lb-in 

2.  weight  (penetration  force)  up  to  15  lb 

3.  rotational  velocity  up  to  2500  rpm 

4 .  si nkage  up  to  4  in 

The  equipment  consisted  of  the  following  items: 

1.  cone  drive  and  test  stand 

2.  control  system  for  cone  drive 

3.  data  recording  system 

Sinkage  is  measured  by  means  of  a  10-turn  potentiometer  through  a 
special  non-slip  wire  running  between  a  reference  point  on  the  base  plate 
of  the  cone  drive  to  the  sample  platform.  It  is  zeroed  for  every  test  by 
using  the  position  crank.  Zero  sinkage  was  taken  when  the  cone  point 
touched  the  soil  surface.  The  velocity  was  measured  by  a  tach-genera tor 
attached  to  the  base  plate. 

All  measurements  were  recorded  on  a  Sanborn  150  strip  chart  recorder. 
The  drive  system  consists  of  an  EC  M0T0MATIC  DC  servo  motor  generator 
with  a  control  panel.  The  capacity  of  the  motor  is  80  oz  in.  In  the 
beginning,  it  was  felt  that  torque  could  be  measured  through  the  control 
system  of  the  motor,  but  this  was  abandoned.  Also  the  original  framework 
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allowed  too  many  vibrations  and  innacuracies  in  measurement  of  the  sinkage  H 
which  is  extremely  critical,  si.icc  H  appears  in  the  third  and  second 
power  in  the  calculations  of  the  shear  stress  ~  and  the  normal  stress  o 
(see  Reference  I  and  Equations  (I),  page  ). 

SOIL  SAMPLES 

The  soil  was  initially  tested  in  round  containers  which  measured  6"  x  6'" 
(large  size  coffee  cans).  The  sand  was  processed  after  each  penetration- 
rotation  test  by  thoroughly  mixing  and  then  shaking  it  a  certain  number  of 
times.  The  grain  size  distribution  is  given  in  Fiqure  k. 

For  the  wet  samples,  the  sand  was  mixed  and  shaken  and  then  allowed  to 
reach  uniform  conditions  by  aging  between  16  to  20  hours.  Samples  were 
taken  before  and  after  every  test  and  moisture  conten*  determined  by  weighing 
the  sample  before  and  after  oven  drying. 

The  shaking  equipment  consisted  of  an  1800  rpm  electric  motor  mounted 
on  a  wooden  plate.  The  motor  shaft  had  a  small  eccentric  weight  attached 
to  it.  The  vibration  thus  generated  did  an  excellent  job  of  compacting 
the  sample.  By  timing  the  duration  of  the  vibration,  different  densities 
could  be  obtained  which  were  very  consistent. 

RESULTS  AND  DISCUSSION  OF  PRELIMINARY  TESTS 

Two  scries  of  preliminary  tests  were  run  during  which  several  improvements 
were  made  to  method  and  equipment. 

Various  moisture  contents  and  densities  were  tested. 

All  tests  were  performed  at  the  basic  standard  velocity  Vq  =  20  rpm. 
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FIGURE  4.  PARTICLE  SIZE  DISTRIBUTION 


The  resulting  data  showed  a  large  amount  of  scatter.  It  was  concluded 
that  the  scatter  in  the  data  was  mainly  caused  by  the  interference  between 
the  cones  and  the  sides  and  bottom  of  the  container. 

Subsequent  tests  were  therefore  carried  out  in  buckets  measuring  10" 
by  10"  wnich  were  filled  with  9"  of  sand.  The  data  thus  obtained  showed 
considerably  less  scatter  for  both  cones. 
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DETERMINATION  OF  STANDARD  a  AND  p 

o  o 

The  standard  values  of  ry  and  p  were  determined  at  a  basic 

o  o 

standard  velocity  of  rotation  ,  Vq  •  This  velocity  was  selected  on  the 
basis  of  the  lowest  practical  and  stable  velocity  the  test  set-up  could 
maintair.  It  was  found  to  be  20  rpm. 

Plots  of  the  t  versus  a  obtained  from  the  experiments  are 
presented  in  Figures  5  and  6- 


where  T  =  the  measured  torque  applied  to  the  cone 
W  =  the  measured  vertical  load  on  the  cone 
H  =  the  measured  sinkage  of  the  cone  tip  below  the  surface 
Cj  iC^  =  constants  dependent  upon  the  apex  angle  of  the  cone ^ ^ 


In  order  to  understand  the  causes  of  the  data  scatter  exhibited  in 
Figures  5  and  6  analysis  was  made  using  a  stepwise  multivariate  regression 
technique  employing  moving  averages  to  fit  a  curve  to  the  experimental  values- 
To  do  this,  averages  of  experimental  values  of  T  and  H  were  plotted  against 
average  values  of  W.  Smooth  curves  drawn  through  these  data  were  then 
used  to  obtain  one  value  of  T  (called  the  curve  value  ,  T_)  and  one  value 
of  H  (Hc)  for  each  value  of  W.  The  ratios 

rt  =  t"  and  *h  =  if 

c  c 
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70  30  ^  A  m  =  25.5 
70305  B  m  =  15.5 


70  30  4  C  m=  25  7 
70  305  D  me  13.3 
70  30  6  D  m  =  I  10 
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were  then  computed  and  plotted  against  each  other- 

Figure  7  shows  three  possible  results  of  such  an  analysis.  In 
Figure  7a  the  conclusion  would  be  that  >s  independent  of  R^ 

i  *  e  -  ,  the  experimental  values  of  T  scatter  about  their  expected  values 
independent  of  the  scatter  of  H  about  H  •  If  the  data  group  as  shown 
in  figure  7b,  it  would  indicate  that,  at  values  of  H  below  the  expected 
value  Hc  ,  the  T  values  are  higher  than  expected,  thus  indicating  an 
inverse  relationship.  The  relationship  depicted  in  Figure  7c  shows  the 
inverse  of  that  shown  in  Figure  7b.  The  scatter  along  the  band,  in  the 
direction  of  the  arrows,  indicates  that  some  other  parameter  influences 
the  phenomenon.  In  most  cases  this  would  be  some  uncontrolled  or  badly 
controlled  test  condition. 

Since  in  our  case  the  variables  which  are  included  in  the  analysis 
are  W  ,  y  and  m  ,  i t  i s  to  be  expected  that  the  data  would  group 
themselves  in  mo i s tu re-dens i ty  r e 1  a  ted  clusters  on  the  R^  and  R^ 
plots.  Similarly,  two  different  cone  shapes  should  show  the  angle  effect 
if  it  exists.  In  Figures  8  and  9  two  test  series  were  carried  out  for  ouch 
cone.  The  loose  soil  was  prepared  by  stirring  and  loosening  air  drv  sand 
in  the  large  sample  can.  ihc  dense  sample  was  obtained  by  using  the 
vibrator  for  25  seconds  on  the  same  air  dry  sand. 

The  density  in  each  case  as  measured  by  a  cone  penetrometer  over 
a  3"  depth  was : 

Loose  :  G  =  4  ps  i  /  i  n 
Dense  :  G  =  I  5  ps i / i n 
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From  the  figures,  it  is  clear  that  scatter  in  the  direction 

is  due  to  density.  For  both  cones  there  is  a  distinct  polarization.  It 
was  to  be  expected  that  the  experimental  sinkage  value  would  be  smaller 
than  the  overall  expected  sinkage  value  for  the  denser  material,  since 
the  overall  value  includes  both  sets  of  data- 

The  scatter  in  the  R.^.  direction  is  thus  a  result  of  other 
influences  on  the  T-H  relation,  since  moisture  is  a  constant  for  both 
the  loose  and  dense  conditions-  It  can  be  seen  that  w i thin  each  density, 
the  torque  T  tends  to  increase  with  respect  to  the  expected  ',rque 
value  ,  ,  when  the  sinkage  H  increased  with  respect  to  the  expected 

sinkage  value  ,  ,  which  is  logical.  The  word  "expected"  should  be 

interpreted  here  in  the  following  manner-  Since  the  values  of  T^  and 
were  calculated  for  the  same  load  W  it  is  o  be  expected  that  the 
actual  data  points  will  spread  themselves  with  some  amount  of  scatter 
about  this  point-  Since  this  point  is  not  really  an  average  in  the  usual 
sense,  it  was  felt  appropriate  to  call  this  central  value  the  expected  value- 

In  Figures  10  and  11  the  versus  R^  plots  for  a  series  of  tests 

at  low  but  varying  moist  ire  content  are  given.  |-  these  tests  attempts 

were  made  to  keep  the  density  constant,  but  it  can  be  seen,  especially 

in  the  60°  cone  test,  that  in  some  of  the  tests,  notably  /0h08D,  the 

density  control  was  poor,  thus  causing  a  more  pronounced  scatter  in  the 

R  direction.  The  30°  cone  data  (Figure  9)  show  that  the  scatter  about 
h 

the  expected  value  of  T(R  =  1-0)  was  independent  of  the  scatter  about 
the  expected  value  of  H(Rh  =  1.0).  Tire  graph  also  sfiows  that  for  the 
depicted  moisture  range  no  conclusions  can  be  drawn  as  to  the  effect  of 
moisture.  This  is  no  doubt  partly  due  to  the  interaction  of  moisture  rod  density. 
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It  was  however  felt  that  the  other  reason  for  this  scatter  was  the 
interference  between  the  container  and  the  cone.  This  interference  would 
also  show  up  differently  for  the  30°  and  60°  cones. 

All  subsequent  tests,  embodying  the  main  part  of  the  investigation 
and  described  in  the  next  chapter,  were  performed  in  enlarged  containers- 
The  resulting  analysis  of  the  data  shovd  a  :harp  reduction  in  the  scatter- 

From  Figures  5  and  6,  the  conclusion  can  be  drawn  that  for  most  of 
the  range  of  values,  the  ielationship  between  t  and  o  is  linear- 

A  regression  analysis  of  the  data  which  was  run  on  the  computer  gave 
the  following  results: 

1 ■  60°  Cone 

Equation  t  =  .001  +  .297  o 

F-Value  =  270 

Degrees  of  Freedom  =86,  1 
Standard  Error  of  Estimate  =  -07 
Correlation  Coefficient  =  -87 

2-  30°  Cone 

Equation  t  =  -  .008  +  .277  cr 

F-Value  =  100 

Degrees  of  Freedom  =  78,  I 
Standard  Error  of  Estimate  =  .100 
Correlation  Coefficient  =  -75 

The  above  data  show  that: 

a)  The  30°  Cone  had  more  scatter  than  the  60°  Cone. 
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EFFECTS  OF  VELOCITY 

Soil  Cond i t i ons 

Tlie  sand  was  kept  at  two  moisture  contents:  air  dry  and  saturated 
(277.  by  weight).  The  density  of  the  air  dry  sand  was  varied  with  the 
vibrator  described  in  Section  3-  The  density  levels  were  defined  as 
Loose,  Medium  and  Dense,  which  conformed  to  volume  weight  ,  y  ,  of  92, 

100  and  1 09  lb/ft^,  respectively. 

The  wet  sand  had  a  constant  density  which  was  determined  by  the 
saturation  moisture  content.  To  attain  greater  uniformity  ,  these 
samples  were  also  vibrated. 

Vc loc i ty  Cond i t i ons 

The  velocities  at  which  the  tests  were  run  ranged  from  about  20  rpm 
to  2000  rpm.  The  velocity  was  set  on  the  DC  control  panel  and  was 
measured  continuously  by  means  of  a  tachometer-generator.  Thus  variations 
of  the  velocity  could  be  observed  during  the  test. 

TEST  METHOD  -  VELOCITY  EFFECTS 

Two  methods  of  testing  were  used  and  the  results  were  analyzed  with 
the  aim  of  determining  whether  different  results  were  obtained- 

Me  thod  I  :  The  soil  was  processed  to  the  desired  density  and  was 
reprocessed  for  each  new  velocity. 

Me  thod  2 :  The  soil  was  not  reprocessed  for  each  velocity.  From 
the  viewpoint  of  ease  of  testing,  the  second  method  is  preferable  since 
the  sample  is  not  taken  out  of  the  tester  and  the  cone  stays  in  place 
over  the  whole  range  of  velocities-  On  the  other  hand,  there  is  no 
knowledge  about  any  change  of  density  during  testing  and  because  of  the 
test  itself.  Only  the  initial  density  is  known. 
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The  first  method  requires  more  time  since,  for  each  test,  the 
sample  has  to  be  reprocessed  or  a  new  sample  has  to  be  inserted-  Some 
amount  of  scatter  will  be  introduced  due  to  the  variation  in  density  between 
the  samples  for  each  test,  as  was  brought  out  also  in  preceding  sections. 

RESULTS  AND  DISCUSSION  OF  VELOCITY  TESTS 
Dry  Sand 

Results  for  the  tests  are  shown  in  Figures  12  through  I/-  They  include 
data  obtained  by  both  test  methods  described  in  the  preceding  section. 

Tli  i  s  inclusion  was  judged  justified  after  analysis  did  not  show  significant 
difference.  Further  proof  was  found  after  the  data  were  plotted,  since 
the  total  scatter  appears  to  be  very  limited  in  magnitude. 

The  first  conclusion  from  the  results  is  that  for  a  steel  cone-sand 
system,  the  relationship  between  i  and  a  is  independent  of  the  velocity, 
since  all  points  fall  along  the  same  curve  and  no  grouping  along  the 
curve  is  evident-  For  the  dry  soils  at  all  three  densities,  and  for 
both  cones,  the  relationship  can  be  expressed  as 

T  -  C,  +  C2  rr 

with  C|  ,  as  expected,  having  a  very  low  value-  The  values  of  the 
constants  for  the  dry  soil,  are  given  in  table  1- 

From  Table  I,  the  following  can  be  concluded: 

a.  The  value  of  C^  for  the  30°  cone  is  very  sm  all-  The 
standard  error  of  the  estimate  is  large  compared  to  Cj  and 
includes  the  zero. 
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air  dry  sand 


AiR  DRY  SAND 


AIR  DRY  SAND 


S  VS  NORMAL  STRESS  RELATIONSHIP  AT  VARYING  VELOC  !  T  I  E  S  -  30°  CONE 


VALUES  F POM  LINEAR  REGRESSION  ANALYSIS  FOR  VELOCITY  TESTS  -  DRY  SAND 
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b-  The  value  of  for  the  60°  cone  is  small  but  appears 

to  be  significantly  different  from  zero. 

c.  The  differences  between  the  C  values  from  within  each 

cone  do  not  give  a  clue  with  respect  to  their  trend-  No  conclusion 

can  be  made  regarding  the  variation  of  C2  with  density. 

d  .  The  difference  between  the  values  for  each  cone  cannot 

be  determined  although  the  average  value  for  the  60°  cone  appears 
to  be  lower  than  the  average  value  for  the  30°  cone- 

Both  figures  and  table  show  that  the  scatter  for  these  series  was 
much  less  than  for  the  preliminary  tests-  This  is  undoubtedly  due  to  the 
experience  gained  in  controlling  test  conditions. 

For  a  dry,  sandy  soil  no  cohesive  or  adhesive  properties,  as  defined 
in  the  usual  sense,  would  be  expected  to  exist.  Thus  the  values  of  Cj 
for  the  30°  cone  would  be  in  accord  with  this  surmise.  It  does  appear, 
however,  that  due  to  the  particular  test  conditions,  some  pseudo-coadhes i ve 
effect  exists  in  the  case  of  the  60°  cone,  since  the  values  of  Cj  here 
appear  to  be  different  from  zero.  Th i s  does  not  necessarily  mean  that 
foi  a  larger  apex  angle,  the  Cj  value  would  increase-  It  is  conceivable 
that  a  maximum  exists  somewhere  between  the  3^°  cone  and  a  flat  plate 
(a  cone  apex  angle  of  180°). 

The  value  of  C ^  is  dependent  by  definition  on  the  material  properties 
of  the  interface  and  thus  it  would  be  expected  that  they  would  be  identical 
for  both  cones.  However,  since  no  particular  care  was  taken  in  the 
polishing  of  the  cones,  some  differences  in  C?  values  may  be  due  to 
surface  roughness  variations- 
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C a  1  culatmn  of  a  gnd 
The  values  of  v 


u-Va lues  and  Their  Velocity  Dependence 

1 

and  u.  were  defined  in  a  preceding  paper 


as : 


_ C1  (2) 

2/3  -  cot  g  (0/2) 

C9  (3) 

n+1  -  c.  cot  g  (0/2) 
n+3  2 


,„  Table  II,  the  values  a!  ,  a  d  u  are  calculated  free,  the  C,  and  C, 
values  of  Table  I.  It  is  evident  that  the  assumptions  made  in  developing 
Equations  2  and  3  are  not  borne  out  by  the  est  results-  There  is  no 
rational  explanation  for  negative  values  for  either  u  or  o  as  obtained 
for  the  30°  cone.  The  principal  reason  for  this  appears  to  be  the 
sensitivity  of  the  denominator  of  both  equations  to  the  magnitude  of  C,  - 

Saturated  Sand 

The  data  for  the  tests  with  saturated  sand  are  shown  in  Figures 
and  19.  u  is  evident  that  a  linear  relation  with  a  positive  intercept 


does  not  exist. 

A  Regression  Analysis  to  fit  a  quadratic  to  the  data  gave  the 

fol lowing  resul ts : 

|  .  30°  Cone 

2 

T  -  -  .033  +  • 2 1  So  +  -003a 

Correlation  Coefficient  -  *98 
Standard  error  of  estimate  =  -055 
F-va 1 ue  =  328 

Degrees  of  Freedom  =  26,2 
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Table  I  I 

CALCULATED  a-  AND  M--VALUES  (n=l) 


Test 

Ang  le 

Cond i t i ons 

Dens i ty 

C1 

lb/ In2 

C2 

a 

lb/ i n^ 

u 

o 

o 

c~v 

Loose 

.01  1 

.22 

negot i ve 

nega  t i ve 

Med i um 

.005 

.24 

I  1 

1  1 

Dense 

.007 

•  174 

1  1 

1.74 

60° 

Loose 

.03! 

.170 

.083 

.373 

Med i um 

.045 

.  132 

.  102 

.253 

Dense 

.036 

.  150 

.088 

.305 
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2 .  60°  cone 

t  =  .021  +  ,091o  +  .09 o2 
Standard  error  of  estimate  =  .011 
Correlation  coefficient  =  .99 
F-value  =  1017 
Degrees  of  freedom  =  21,2 

the  intercepts  again  are  small  and  the  standard  error  of  the  estimate  is 
rather  large  by  comparison  although  not  as  pronounced  for  the  60°  cone. 

In  general,  the  data  from  the  60°  cone  tended  to  show  less  scatter. 
During  the  tests,  it  was  observed  that  the  60°  cone  tests  were  easier  to 
carry  out  because  of  smaller  penetration.  From  these  tests  no  definite 
conclusions  can  be  made  regarding  the  character  and  interpretation  of  the 
three  constants.  More  data  ove'  a  wider  range  of  moisture  contents  are 
requ i red . 

CONCLUSIONS 

From  the  tests  with  wet  and  dry  sand  the  following  conclusions  can 
be  formulated: 

1)  There  appears  to  be  a  linear  relationship  between  the  shear 
stress  on  a  normal  stress  for  the  cone  in  sand. 

2)  The  assumptions  made  for  the  development  of  equations  2  and  3  are 
not  borne  out  by  the  test  results;  this  is  partly  due  to  the  sensitivity 
of  those  equations  to  the  slope  factor  C ^  of  the  test  data. 

3)  The  relationship  between  t  and  a  does  not  appear  to  be  influenced 
by  velocity  for  a  sand  in  the  dry  or  saturated  condition. 
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RECOMMENDATIONS 

It  is  recommended  that  funds  be  made  available  to  carry  out  further 
investigations  of  the  following  points: 

1.  A  wider  range  of  moisture  content  for  a  sand  in  order  to 
define  more  clearly  the  parameters  obtained. 

2.  Interface  parameters  and  their  velocity  dependence  for  cohesive 
and  mixed  soils  at  various  moisture  contents. 

3*  Interface  parameters  and  their  velocity  dependence  for 
coated  cones  in  various  soils  at  various  moisture  contents- 

4.  Performance  of  tests  with  wheels  and  reduction  of  data  using 
dimensional  analysis  methods  in  order  to  incorporate  the  interface  parameters 
into  a  prediction  equations  system. 
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Appendix 

1.  SORTING  PROGRAM  (Only  for  Themis) 


D  X  ME-Jb  I  O'l  w(50> , IVS(2000> , VEL(2000) , 

2TA'J(?0?I0>  ,SIGMA(20C0> 

INTEGER  TE^TNQ 

CALL  iriLE<20, 'CONEV*  ) 

call  ur i le; < 21  *  'Cdut'  ) 

1*0 

11=1 

HE40(2^,?A)  K A  »  KP . KM , KV 

2^  KORMATMI) 

32  KLAU  (20, 20)  TESTNO, DIG, ID.SOIL, THETA, DENS. WATER 
IF  (EOKP  (2W)  )  GO  TO  30 

20  FORMAT ( I?, Al, II, A4.F3.0, A3.F) 
if  (Test mo, to, ^(0)  go  to  30 
IF  (THFTA.EU.30, >  K«KA 

IF  (THETA,  FCI,  60,  )  K  a  2  *  K  A 
IF  (OEMS, EG. 'DNS' )  K*K+KD 
IF  (OEMS,  ED.  'MFD'  )  K=K+2#KD 
IF  (OEMS.EO. »LSE' )  K*K+3»KD 
IF (WATEH,En,0,0)  KrK+KM 

IF  (WATER. GT. fa, 0.ANU, WATER, LE,  5,0)  K=K+2*KM 
IF  (wATER.GT.5.0. AND, WATER, LE. 10, >  K*K+3*KM 
IF  (WATER. CT. 10, .AND, WATER, LE, 15, )  K»K+4«KM 
IF  (WATER. GT. 15, .AND, WATER, LE, 20, )  K»K+5«KM 
IF  (WATER, GT. 20, .AND, WATER, LE, 25,  )  K«K+6*KM 
IF  (WATER. GT, 25, .AND, WATER, LE,30.  )  K»K  +  7*KM 
IF  (WATER, GT.30|.ANU,WATER,LE, 35,  )  K*K+8*KM 
IF  (WATER, GT.35, .AND, WATER, LE, 40.  )  K*K+9*KM 
TAMG=5INP(THETA/2, )/C0SD<THETA/2, ) 

StCA  =  l./C0SD(THETA/2,  ) 

31  HLAD(27',21)  T  ,  W  |  V  ,  H 

21  FORMAT ( 4  F  ) 

IF  (T.FQ.0,  )  GO  TO  33 
1  =  1+1 

TAU{ I )=T/(3,142*(TANG**2>*SECA+(H**3)  ) 

S I GMA  C I )=W/(3,142*TANG##2)*(1,/(H**2) ) 

IVS( I )  =  I  NT ( V ) *KV*K  +  I  I 
VEL(  I)=V 
GO  TO  31 

33  N (  I  I  )  =  I 
I  I  8  I  I +1 
GQ  TO  32 
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30  NDAT=I 

1)0  41  I»1,NDAT-1 
□0  41  Jal.NDAT 

IF  ( I V9( I ) ,GT,  I VS( J)  )  GO  TO  41 

ITlamU  ) 

T2=TAU< I  ) 

TJ  =  SIGHA(  I  > 

T4  =  V£L(  1  > 

IVS{ I )*IVS(J) 

T A U ( I ) =T  AU ( J  ) 

S  J  GM A ( I ) =S I GMA ( J ) 

VtL ( I ) s VFL  <  J ) 

I  VS  < J>  =  IT1 
TAU(  J)=>T? 

S I GM A ( J ) =T3 
VEL( J)=T4 
41  CONTINUE 

WHITP  (21*22) 

E.MU1  =  TAU(1>/SJGMA<1> 

UO  42  Isl,MDAT 


EMUaTAUd  >/SIGMA<  I  ) 

KmUsEMU/EMUI 
HVeL=VFL( I )/VElCl) 

42  WHITE(2i;23)  I  VS ( I) * T AU ( I ) * S I GM A ( I ) . EMU  * RMU , R VEL 
23  F0RMAT(3X* I12*5F9,4) 

22  FORMAT (1H1* ’ HEA0ING  COMES  HERE'*5(/)) 

CALL  E.V  I  T 
END 


The  program  will  sort  according  to  the  moisture  content,  density,  cone 


angle,  and  velocity.  The  sorting  is  done  according  to  a  code,  defined  by 
the  input  constants  KA,  KD,  KM  and  KV.  The  magnitude  of  these  constants 
determines  the  magnitude  and  location  of  the  variable. 


For  instance: 

The  number  1011012018  Is  obtained  by  setting 

KA  =  1.000,000,000 

KD  =  1,000,100 

KM  =  10,000,000 

KV  =  100 
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and  signifies  that  this  test  is  for  a  30°  cone,  0.01  moisture,  dense  sand, 
velocity  of  120  rpm  and  belongs  to  test  set  number  18. 


2.  REGRESSION  PROGRAM 


DIMENSION  X  <  4fcJPl )  ,TAUU00>  *SIG<100)  ,STD<5>  .XBAr<5>  ,D(10) 

♦  » ISAVE(5) »RX(16) * RY(4)»LV(4) ,MV(4)i B(4) tSB(4)» T(4) 

♦  i  ANS  ( 1 0  )  ,R(10> 

ID  1=1 

TYPE  3 

3  FORMAT!*  ENTER  INPUT  BLOCK'/) 

2D  ACCENT  1,TAU( I ) ,SIG(  I  ) 

IP (TAU( I ) .EU.99JG0  TO  21 
1  =  1*1 
GO  TO  20 
21  N=  1  -1 

DO  23  I  = 1 »  M 
X( I  )  *S  I G (  I  ) 

X( I*N)*SIG(  I  )#*2 
X( I*2*H)eSIG(  I  )**3 
23  X( I *3#N )  =  T AU (  I  ) 

M  s  4 
10  =  1 

CALL  C3RRE(N,M»I0iX#XBAR»STDjRX,R#D>BiT,i 
TYPE  2 i ( V  < I ) ,X( l*N> . X  C I*2*N) ,  X< I*3*N) , I?1,N) 

DO  3P  K  =  1  *  3 
DO  25  Jsl.K 
25  I S A  VE ( J  )  = J 
K  P  1  3  K  ♦  1 

CALL  0RDrR ( KP1 »  R  »  4  *  K  * l S A VE , RX , R Y  ) 

CALL  MNV(RX#K*DET#LV*MV) 

CALL  M'JLTR(N,K,XPAR.STD,D,RX*RY,  ISAVE.B.SB.  T,  ANS) 

TYPE  4,  (ANSI  I)  ,  1=1,10)  ,  (B( I ) ,SB( I ) ,T(I ) , Jal,K) 

3D  CONTINUE 
GU  TO  10 

1  FORMAT (2r) 

2  FORMAT  <34 { 4F/  )  ) 

4  FORMAT (4F/4F/2E/3(3F/J  ) 

LNO 

subroutine  data 
return 

ENO 


This  program  is  set  up  to  be  used  on  the 
computer.  Data  are  entered  in  pairs  of  t,o 


teletype  with  a  PDP-10 
on  the  teletype. 
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The  subroutines  Corre,  Order,  Minv  and  Multr,  are  standard  SSP 
subroutines  -  see  IBM  application  Program  H20-0205-3* 

The  output  of  this  program  is  as  follows: 

2  3 

1.  Input  data  are  repeated  and  a  and  a  values  are  printed  out. 

2.  ANS(I)  1=1,  10 

ANS  I  =  Intercept 

2  =  Multiple  Correlation  Coefficient 

3  =  Standard  Error  of  Estimate 

4  =  SSAR  (Regression) 

5  =  DF  Associated  with  SSAR 

6  =  Mean  Square  of  SSAR 

7  =  SSRD  (Deviations  from  regression) 

8  =  DF  Associated  with  SSDR 

9  =  Mean  Square  of  SSDR 
1 0  =  F-Val ue 

3-  B(l),  SB ( I ) ,  T ( 1 )  I  =  1,  K 

B(l)  =  first  column  at  the  regression  constants 

SB ( I )  =  second  column  -  the  standard  duration  associated 
with  B(l) 

T(l)  =  third  column  -  the  t-value  associated  with  B  ( I ) 
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